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ARE THERE ENOUGH IONIZING PHOTONS TO REIONIZE THE UNIVERSE BY Z « 6? 
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Draft version February 2, 2008 

ABSTRACT 

An estimate for the number of ionizing photons per baryon as a function of redshift is computed based on the 
plausible extrapolation of the observed galaxy UV luminosity function and the latest results on the properties 
of the escape fraction of ionizing radiation. It is found that, if the escape fraction for low mass galaxies 
(Mtot S 10 n M Q ) is assumed to be negligibly small, as indicated by numerical simulations, then there are not 
enough ionizing photons to reionize the universe by z = 6 for the cosmology favored by the WMAP 3rd year 
results, while the WMAP 1st year cosmology is marginally consistent with the reionization requirement. The 
escape fraction as a function of galaxy mass would have to be constant to within a factor of two for the whole 
mass range of galaxies for reionization to be possible within the WMAP 3rd year cosmology. 

Subject headings: cosmology: theory - cosmology: large-scale structure of universe - galaxies: formation - 
galaxies: intergalactic medium 



1. INTRODUCTION 

Studies of cosmic reionization - especially theoretical ones 
- have never been considered as a "photon-starved" field. The- 
orists always felt free to adjust the emissivity of the ionizing 
sources, usually quantified by the escape fraction of ionizing 
radiation, to adjust the reionization redshift to their choosing. 

This approach was, indeed, justified in the earlier stud- 
ies, since until recently limited knowledge existed on the 
reasonable values for the escape fract ion at high redshifts . 
However, the latest observational dGiallongo et al] l2002t 
Fernandez-So to et al l l200l IShaplev et all 120061: IChen et al I 



2007) and numerical (|kazoumov & Sommer-Larsenr t2006, 
20071: [Gnedin et al. 2007) studies finally begin to converge 
on the values and evolution of the escape fraction of ioniz- 
ing radiation and on the relative escape fraction between the 
far UV and Lyman limit. Three properties of the escape frac- 
tion are particularly important for reionization studies: (1) 
the value of escape fraction is small (a few percent at most), 
which is an order of magnitude smaller than is assumed in 
some reionization modeling, (2) it is weakly dependent on the 
galaxy mass or star formation rate for large galaxies, and (3) 
it drops prodigiously for dwarf galaxies. The first two prop- 
erties are reproduced in all recent studies, both observational 
and theoretical, and, therefore, are rather robust. The last fea- 
tu re of the escape fra ction has only been seen in simulations 
of iGne din et al. (2007 |) and is indicated by measurements of 
iFernandez-Soto et al.l (12003I) 4 . because other simulations and 
observational studies do not yet have either numerical resolu- 
tion or sensitivity to resolve dwarf galaxies. 

Another important observational advance that places the 
study of reionization on a much more quantitative footing is 
the observational determination of the galaxy UV luminosity 
function down to well below at z > 6. Since it is not possi- 
ble to give a comprehensive review of all observations in this 
Letter due to space limitations, I refer the reader to the recent 
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work by iBouwens et al.l (120071) . who give a detailed review 
of the current status of existing observational data. While the 
data for the galaxy luminosity function during the reionization 
era (z > 6) are still sparse, the plausible extrapolation of the 
observed z ~ 6 luminosity functions to higher redshifts can be 
used to predict the global production of the ionizing radiation 
to at least within a factor of 2 to 3, i.e. more than an order of 
magnitude improvement over the previous, purely theoretical, 
assumptions. 

Of course, the most accurate models are only possible with 
the large-box and high-resolution cosmological simulations, 
which model in detail the emission of ionizing photons in high 
redshift galaxies and quasars, the propagation of ionizing ra- 
diation in the expanding universe, and absorption of that ra- 
diation at cosmic ionization fronts and Lyman limit systems. 
But even the simplest balance of the available ionizing pho- 
tons and the number of atoms that need to be ionized before 
the end of reionization at z ~ 6 - as required by the observed 
transmitted flux in the spectra of high redshift qu asars discov- 
ered by the SDSS collaboration dFan et a l. 2006) - is a useful 
exercise after the recent improvements in our understanding 
of the sources of reionization. 

This is the subject of this Letter. 

2. RESULTS 

In order to compute the total number of ionizing photons 
available for reionizing the universe at any given redshift, the 
observed luminosity functions need to be extrapolated to ear- 
lier redshifts. Such extrapolation is, of course, not unique. 
However, since the mass function of dark mater halos can be 
computed sufficiently precisely in a given cosmology at any 
redshift, the extrapolation of the luminosity function to z > 6 
can be made reliably if the relationship between the galaxy 
luminosity and the mass of its dark matter halo can be estab- 
lished. 

While such a relationship is unlikely to be a simple 
function, models that assume a one-to-one correspondence 
between the galaxy luminosity and the halo mass pro- 
vide remarkably go od fit to a variety of observational tests 
dConrov et alj|2006l) . Thus, as a simple and crude approxima- 
tion, it is instructive to assume such a relationship between 
the galaxy UV luminosity Lj/v for the high redshift galaxies 
as well. 
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FIG. 1 . — The cumulative mass function of dark matter halos (top panel), 
the cumulative UV luminosity function of galaxies (right panel), and the re- 
lationship between the UV absolute magnitude and the halo mass (central 
panel) obtained by matching the two. Four different lines correspond to 
Z = 3.8 (solid lines), z = 5.0 (long-dashed lines), z = 5.9 (short-dashed lines) 
and z = 7.4 (dotted lines). The obser vational data are repre sented by their 
respective Schechter function fits from Bouwens et al. (2007). 




FIG. 2. — The relationship between the UV absolute magnitude and the 
halo mass from Fig.fT] shifted horizontally and vertically using the correction 
factor C = (1 +z)/7 (line styles are as in Fig.ff}. Notice that the shift results in 
a good agreement between different redshifts for all galaxies more massive 
than about 3 X 1O 1O M0. The blue/black solid line marks the average mass- 
to-light relationship used in the rest of this paper. 

In a given cosmology, this relationship between the total 
mass of a dark matter halo Mtot and the luminosity of a 
galaxy hosted in that halo Luv can be obtained by matching 
the cumulative mass and luminosity functions, 



«(> Luv) = «(> Mtot), 



(1) 



as is shown in Figure [T] for the 4 values of redshift for which 
iBouwens et al.l d2007l) give the parameters of the Schechter 
function fits, assuming the best-fit cosmology for the combi- 
nation of WMAP 3rd year data and Large Re d Galaxies part of 
SDSS survey (WMAP3: ISpergel et al.ll2006l) 5 . Unfortunately, 
the derived relation between the UV luminosity (as expressed 
by the AB absolute magnitude at 1600, .A/fi6oo) and the total 
halo mass is redshift-dependent, and so cannot be easily ex- 
trapolated to higher redshifts. However, a simple correction 
of the UV luminosity by a factor of C and the total mass by a 



C(Z): 



1+z 
7 



eliminates most of the redshift dependence for halos more 
massive than about 3 x 10 10 M Q , as is shown in Figure[2] As 
I discuss below, these low mass halos contribute almost noth- 
ing to the ionizing photon budget, and are so unimportant for 
the purpose of this paper. Throughout the rest of this paper, I 
use the average relation marked by a blue/black line in Fig. [2] 
The relation between the UV and ionizing luminosi- 
ties is quantified by the relative escape fraction, / esC| rei = 
/esc(912)// esc (1600) and the value of the intrinsic ratio of stel- 
lar luminosities at these two wavelength, r; nt = (Li60o/L9i2)int- 
F or the escape fractio n at the Lyman limit I adopt the results 
o flGnedin et al.l(l2007l) . who found that in high resolution sim- 
ulation of galaxies with radiative transfer escape fractions for 
larger galaxies are of the order of a few percent, consistent 
with observational determinations, but little (if any) radiation 
escapes from small galaxies. Thus, for the relative escape 
fraction as a function of galaxy mass I adopt the following 
form: 



/esc,rel(M)«0.15 



1 

^min 



if Mjot > 5 x 1O 1O M , 
otherwise. 



This form is consistent with o bservational measu rements of 
the relative escape fraction (c.f. Shaplev et al. 2006) for mas- 
sive galaxies. The drop in the escape fraction is also indicated 
in observations of Fernan dez-Soto et"ai1 d2003l) . but the char- 
acteristic transition occurs at a factor of 10 higher star forma- 
tion rate, which would correspond to a higher characteristic 
mass. The adopted value of 5 x 10 10 M Q therefore likely bi- 
ases the estimate production of ionizing photons up. Here, 
as well as in the rest of the paper, all uncertain quantities are 
chosen so that to insure that my estimate for the total number 
of ionizing photons is likely to be an overestimate, rather than 
an underestimate. I return to the uncertainty of the main result 
in th e Discussion section . 

In lGnedin et al.l (120071) simulations, the relative escape frac- 
tion of low mass galaxies (Mtot < 5 x 1O 1O M ) never ex- 
ceeds about 0.01, and is often much lower, so I adopt 5 m ; n = 
0.05 as my fiducial value. I consider the effect of parameter 
s m i„ in the Discussion section. 

Using the mass-to-light matching from Fig. [2] the mass de- 
pendence can also be recast as the luminosity d ependence. 

For the intrinsic break, I adopt a value of dShaplev et al.l 
120061) 

r 'mt ~ 3. 



However, Siana et al] (12007b argue for a larger value for the 
intrinsic break, r mt w 6, and that larger value is also consistent 
with Starburst99 spectral synthesis models ( Leither er et al.l 
1999). Here I again adopt a lower value as a fiducial num- 
ber so that not to underestimate the total number of ionizing 
photons. 

With the above assumption, the total emission rate density 
of the ionizing photons at a given redshift z can now be ex- 
pressed as 



n 7 = / dL\ 



esciel M600 



dn 



1600- 



r m (E) dL 



(2) 



1600 



5 n„ = 0.27, h = 0.71, n s = 0.95, <r 8 = 0.78 



where (E) is the average energy of a photo-ionizing photon 
(which I ta ke to be 22 e V, consistent with typical spectra of 
star-bursts dLeitherer et alj [l999l) ) and dnjdLium is the co- 
moving UV luminosity function obtained from the halo total 
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FIG. 3. — The total ionizing photon-to-baryon ratio as a function of 
redshift for the WMAP3 (red/gray lines) and WMAP1 (blue/black lines) 
cosmologies. In both cases the dashed lines show the contribution from 
galaxies only, while solid lines include the quasar contribution (with 
J min = 0.05). The shaded region is an estimate of the upper 95% CL 
uncertainty (the lower uncertainty is harder to estimate, and so its exact limit 
is not shown). Two black dotted lines mark the two estimates of the quasar 
contribution. The adopted reionization criterion, 1 < N^ij, < 3 at z = 6, is 
shown as a segment with error-bars. 

mass function using the mass-to-light ratio from Fig. [2] The 
total number of ionizing photons per baryon at time t is then 

N l/b {t)= [ ^dt, 
Jo n b 

where «/, m 2.5 x 10~ 7 cm" 3 is the comoving number density 
of baryons. 

An additional complication in this estimate, however, is 
introduced by a possible contribution to the ionizing back- 
ground from high-redshift quasars. Unfortunately, this con- 
tribution cannot be estimated with the same method as the 
contribution from galaxies, because quasars are known to be 
short-lived, and halos of a given total mass may or may not 
host a quasar at any given time. 

Fortunately, the quasar luminosit y function is reaso nably 
well known all the way to z ~ 5 (Hopkins et al. 2007). Us- 
ing the fitting code provided by Hop kins et al.1 (120071) . the 
emission rate density of ionizing radiation from quasars can 
be estimated at a range of redshifts 4 < z < 6. Such an es- 
timate agrees rema rkably well with an earlier estim ate by 
iMadau el : all (1 1999b : for example, iHopkins et ail (120071) fit at 
z = 5 results in log(« 7 ) « 50.5, while lMadau et al.l (1 19991) esti- 
mate at that redshift is log(/i 7 ) « 50.6. At z = 4 both estimates 
give the same value of log( n-y) » 50.8. 

The extrapolation of the Hopki nset al.l d2007) luminosity 
function to z > 5 is, of course, highly uncertain. In order to 
approximately account for possible uncertainties, I consider 
two different extrapo lations: the "lower" one simply uses the 
Hopkins et al. (2007) best-fit model to compute the quasar lu- 
minosity function at any redshift; the "higher" extrapolation 
multiplies the "lower" one by a factor of [(1 +z)/6] 3 (chosen 
somewhat arbitrarily), greatly increasing quasar abundance at 
higher redshifts. 

The resultant photon-to-baryon ratio is shown in Figure 
[3] for two adopted sets of cosmological parameters: the 
WMAP3 cosmology introduced above and the best-fit values 
for the pur e ACDM model fro m the first year WMAP data 
(WMAP1; Sper gel et ai1l2003l) 6 . As can be seen, the quasar 
contribution is smaller than the one from galaxies by at least 

6 fl M = 0.27, h = 0.72, n s = 0.99, <r 8 = 0.90 



an order of magnitude, and so its large uncertainty is not that 
important. 

The uncertainty on the estimates from Fig. |3]is not easy to 
evaluate. The observational 2a limit on the adop t ed value of 
/ esc .re i « 0.15 is about 50% (IShapley et al.ll2006h IChen et alJ 
2007). In addition, the model for the evolution of the ioniz- 
ing luminosity from galaxies presented above is only approx- 
imate. Since the observed luminosity functions are only used 
to compute the total number of ionizing photons (an integral 
quantity), the observational 2a error f rom the uncertainty i n 
the luminosity functions is about 40% (Bouw ens et al. '2007). 

1 increase it to 50% because additional interpolation is in- 
volved, and adding the two uncertainties in quadrature, ob- 
tain about 70% uncertainty in the upper direction. The lower 
uncertainty is much larger, since the selected values for all 
relevant factors are consistently biased toward a larger value 
of A/ 7 //,. For example, the value for the intrinsic break n nt 
may be close to 6 than to 3, resulting in additional factor of 

2 uncertainty in the downward direction. The lower shaded 
regions in Fig. [3] include this factor-of-2 uncertainty, but the 
lower bound of the uncertainty is likely to be even lower, so it 
is not shown in Fig. [3] 

Thus, it is unlikely that curves in Fig. [3] underestimate the 
correct value by more than a factor of two, although they can 
present a substantial overestimate of the correct result. 

3. DISCUSSION 

But what is the reionization criterion in terms of N 1 /i,l 
Obviously, the whole process of reionization cannot be sup- 
planted by a single value for Nj/i,. On the other hand, having 
N-y/b > 0.9 is a necessary condition for reionization (assuming 
that helium is singly ionized and the mass fraction of the neu- 
tral hydrogen at z ~ 6 is small, consistent with the expected 
abundance of Lyman limit systems where most of neutral hy- 
drogen resides after reionization). 

On the other hand, a ratio of the the Hubble time to the aver- 
age recombination time in a particular region of the universe 
is 

^-^I3^y\l + (S))C R , (3) 

where (S) is the average overdensity in that regio n, and C R is 
the r ecombination clumping factor in that region dKohler et al. 
|2003), C R = (R(T)n e n ml )/((R(T))(n e )(n ml )). 

Since th e definition of the es cape fraction of ionizing radi- 
ation from lGnedin et al. (2007) accounts for all local absorp- 
tion, including high density gas inside a galaxy halo, which 
dominates the clumping factor, Cr cannot be large in the gen- 
eral IGM. More than that, the conclusion that reionization is 
complete by z ~ 6 comes primarily from the observations of 
the SDSS quasars dFan et alj [2006, and references therein). 
The transmitted flux at z > 5.5 comes mostly from the centers 
of large voids, where (1 + (S)) f=a 0.1 and Cr < 10. Equation[3] 
then implies that less than 2 ionizing photons per baryon (out- 
side the virial radii of ionizing sources) are needed to satisfy 
the observational requirements. Of course, if the local absorp- 
tion inside the virial radius (characterized by a large clumping 
factor Cr) is included, the required number of ionizing pho- 
tons per baryon will be much higher; but then a correspond- 
ingly larger value for the escape fraction (which excludes lo- 
cal absorption) should be adopted. 

That estimate is also consistent with the conclusion by 
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FIG. 4. — The total ionizing photon-to-baryon ratio as a function of red- 
shift for the WMAP3 cosmology for different assumed values of the param- 
eter s m j n : (dotted line), 0.05 (short-dashed line, the fiducial value from 
Gnedin et al. 12007) simulations), 0.3 (long-dashed line) and 1 (solid line). 

iMiralda-Escudel (120031) . who estimated that 

1 dN ilb < ? 
H dt ~ 

for 6 < z < 9, which translates into iV 7 «, < 2.5 if the contri- 
bution of sources beyond z = 9 is unimportant. 
Thus, a requirement 

K/V 7/i (z = 6)<3 (4) 

appears to be a sensible criterion for the reionization of the 
un iverse by z = 6. The same c ondition has also been obtained 
by iBolton & Haehnell d2007l) from extrapolating the produc- 
tion rate of ionizing photons required to fit the observed evo- 
lution of the mean opacity of the Lyman-alpha forest to z = 6. 

Thus, the WMAP1 cosmology is marginally sufficient to 
satisfy the condition (0|, while the WMAP3 universe is well 

7 But not the top-heavy IMF, since the conclusion presented in this paper 
is independent of the stellar IMF and does not require any assumption about 
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